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The potential between static-light mesons forming a meson-meson or a meson-antimeson system is calculated 
in quenched and unquenched SU(3) gauge theory. We use the Sheikholeslami-Wohlert action and statistical 
estimators of light quark propagators with maximal variance reduction. The dependence of the potentials on the 
light quark spin and isospin and the effect of meson exchange is investigated. Our main motivation is exploration 
of bound states of two mesons and string breaking. The latter also involves the two-quark potential and the 
correlation between two-quark and two-meson states. 



We have been working on lattice calculations of 
multiquark systems in order to understand their 
properties from first principles. Here "multi" 
means that the system can be decomposed into 
more than one colour singlet, the simplest case 
being four quarks. Previously, we have studied 
four quarks in the static approximation in SU(2), 
getting energies for a general set of geometries 
which have been reproduced by a model based on 
two-body potentials and multiquark interaction 
terms [1]. We have also looked at the flux distri- 
bution corresponding to the binding energy [2]. 
These results for the static case are discussed in 
last year's proceedings. 

In the current work only two of the quarks are 
static and SU(3) is used, also unquenched. The 
binding of this system of two heavy-light mesons 
is studied as a function of the heavy quark sep- 
aration. Different cases of spin and isospin of 
the light quarks are measured and the effect from 
meson exchange is extracted from the potentials. 
The model for energies of static quarks is ex- 
tended for this more dynamic system as discussed 
in Ref. [3] and by A.M. Green in these proceed- 
ings. The meson-antimeson case with light quark 
isospin I q = is relevant to string breaking. 

In the future we plan to explore breaking of an 
excited string, which is interesting from the point 
of view of hybrid meson phenomenlogy. 
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1. TWO HEAVY-LIGHT MESONS 

Experimental candidates for bound states of 
four quarks (of which two are antiquarks) lie close 
to meson-antimeson thresholds. In our case e.g. 
a heavy T particle could be a £?*5* system. Sys- 
tems with heavy quarks should be more easily 
bound as the repulsive kinetic energy is smaller 
with the attractive potential being flavour inde- 
pendent. The binding of four-quark systems has 
been studied e.g. in bag, string-flip and deuson 
models, states with two b quarks being stable in 
most models. 

In the lattice calculation we measure two 
diagrams for both meson-meson and meson- 
antimeson cases - see Fig. 1. One of them is the 
unconnected one without light quark interchange 
and the other connected, where the light quarks 
hop from one meson to another. The light quark 
mass we use is approximately that of s. Because 
of the static approximation for the heavy quarks 
their spin and isospin decouples making the pseu- 
doscalar B and the vector B* degenerate, while 
physically they have a 46 MeV (1%) separation. 
We call this degenerate set B. For the two- meson 
system combinations of B and 5* do have differ- 
ent energies in our case. We measure wavefunc- 
tions symmetric under interchange of the mesons 
with the light quark spin and isospin being sin- 
glet or triplet; these then couple to B,B* combi- 
nations [4]. 

We obtain estimates of light quark propaga- 
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Figure 1. Diagrams for BB and BB systems. 

tors using maximally variance reduced pseudo- 
fermionic ensembles [6] , 24 of which are generated 
for each gauge configuration. A second nested 
Monte Carlo calculation is then performed for the 
pseudofermions of each gauge configuration. The 
variance reduction requires "halving" of the lat- 
tice into two separate regions with the propaga- 
tors going from one region to the other. Thus 
the connected diagram for BB and the cross cor- 
relator between QQ and BB need to have one 
of the normally spatial axes as the temporal 
axis, introducing considerable technical compli- 
cations. The parameters of our calculation are 
(3 = 5.2, Csw = 1-76, a « 0.14 fm, M PS /M V = 
0.72 for unquenched [5] and (3 = 5.7, Csw = 
1.57, a « 0.17 fm, M PS /M V = 0.65 for the 
quenched case with a 16 3 x 24 lattice being used 
for both. We have 20 and 54 gauge configura- 
tions for quenched and unquenched respectively; 
with pseudofermionic fields these take some 60 
GB of diskspace. A variational basis of local and 
fuzzed mesonic operators is used, diagonalization 
of which maximizes overlap with the ground state 
of the system. 

We measure the different spin and isospin com- 
ponents as discussed in Ref. [4]. For the BB 
case I q = 1 has only the unconnected diagram, 
whereas I q = has the connected one subtracted. 

2. RESULTS 

The raw correlators show that for the BB sys- 
tem the unconnected diagram is much noisier and 
does not contribute to the binding at R > 1. The 
connected diagram, on the other hand, gives a 
small binding for larger R and also contributes to 
the observables where the spin of the light quarks 
changes. 

At R = the heavy quarks are at the same 



point and the BB case looks like a baryon with 
an antitriplet string. We can compare to previ- 
ously measured [6] energies of the A5, baryons 
for I q , S q = (0, 0) and I q ,S q = (1,1) respectively, 
finding excellent agreement [4] . States with a sex- 
tet string (0,1; 1,0) lie higher. The BB singlet at 
R = looks like a pion, and we find agreement 
with the energy of a pion with non-zero momen- 
tum. 

The meson-meson potentials are shown in Fig. 
2. The I q ,S q = (1,1) case is similar to (0,0) 
but less bound; the level ordering observed at 
R = is retained and the attraction disappears 
earlier. For the (1,0) a remnant of the sextet 
string makes the small- R potential repulsive. The 
(0, 1) at R = is attractive for unquenched and 
repulsive for quenched; this is the only qualitative 
2.5a) difference between the quenched and un- 
quenched results visible in our data. Both (0, 1) 
and (1,0) seem to have attraction at r « 0.3 fm, 
which is a meson exchange effect as opposed to 
the small distance behaviour governed by gluonic 
effects. From a crude two-body Schrodinger ap- 
proach using these potentials we expect binding 
for all of these cases except perhaps (1,0). 

The meson- antimeson potential for I q ,S q = 
(1,0) is shown in Fig. 3, the (1,1) case being sim- 
ilar. For both of these a QQ + qq state with the 
same quantum number is lighter for small R. For 
(1, 0); (1, 1) the relevant energies are V(R)+7r and 
V(R) + p respectively, where V(R) is the static 
QQ potential. An estimate of V(R) + tt is in- 
cluded in Fig. 3. 

The contribution from meson exchange can be 
examined e.g. by looking at the crossed dia- 
gram measured in the quenched approximation 
at fixed T as a function of R. The BB BB 
case should have a contribution from p exchange, 
and we indeed find agreement by using the pre- 
viously measured m p and normalizing by hand. 
For BB* — > B*B we should have a contribution 
from 7r exchange. In this case we can use a re- 
cent determination of the BB*tt coupling [7], the 
experimental decay constant and our m^. In the 
one-7r exchange formula everything is thus known, 
and we find excellent quantitative agreement for 
R > 0.5 fm. This is strong support for deuson 
models [8] in this distance range. 
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Figure 2. BB potentials. Here 2M B Ro = 
4.98(1), 5.83(3) for quenched, unquenched re- 
spectively and aRo = 0.49 fm. 

3. STRING BREAKING 

In our quenched and unquenched calculations 
the ground state BB and QQ potentials cross at 
r w 1.2 fm. We are investigating the break- 
ing of the QQ string by using a variational ap- 
proach similar to that used in Higgs models by 
several groups. The cross correlator between two- 
meson and two-quark states allows us to study 
their mixing also in the quenched theory - in the 
unquenched case additional fermion bubbles in- 
duce corrections. The quenched mixing matrix 
element can then be used to estimate the splitting 
of energy levels at the string breaking point, even 
though no actual splitting occurs with quenching. 
With an unquenched calculation the energy split- 
ting can be studied directly using the full varia- 
tional approach. 

One might think that an excited string would 
break at a smaller distance than the ground state. 
This is not necessarily the case, as e.g. the first 
excited state has J z = 1 with quark separation 
along z and only breaks into mesons BlBu with 
L + V > 0. In general it is an open question if 
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Figure 3. BB potentials. For small R a state with 
a QQ and a pion should be lighter. 

a state with particular quantum numbers has the 
lowest energy at a particular heavy quark separa- 
tion as a) a hybrid QQ meson with excited glue, 
b) a ground state QQ meson and a qq meson or c) 
two heavy-light mesons. These energy levels and 
their mixing can be studied on the lattice with 
our techniques. 
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